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25 ' We present a measurement of the branching fraction of neutral B meson decaying to final states 

containing a Ki meson, i.e. -ftri(1270) and -ftri(1400), and a charged pion. The data, collected 
with the BABAR detector at the Stanford Linear Accelerator Center, represent 454 million BB 

^ ' pairs produced in e"'"e~ annihilation. We measure the branching fraction B{B^ -^ K^it~) = 

H . (31.0 lb 2.7 lb 6.9) X 10^^, where the first error quoted is statistical and the second is systematic. 

In the framework of the i^— matrix formalism used to describe these decays, we also set limits on 
the ratio of the production constants for the 7^1(1270)"'' and -fCi(1400)''' mesons in B^ decays. 
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1 INTRODUCTION 

B meson decays to final states containing an axial-vector meson and a pseudoscalar meson have 
been recently studied both experimentally and theoretically. Branching fractions of the B mesons 
decays to final states containing an ai(1260) or bi meson associated with a pion or a kaon have 
been measured experimentally [1]. Theoretical predictions for the branching fractions of B decay 
modes to final states containing an axial-vector and a pseudoscalar meson have been calculated 
assuming a naive factorization hypothesis [21 13] and QCD factorization [1] . Expected branching 
fractions of these B meson decay modes are of the order of 10~^. 

Recently the BABAR Collaboration has measured CP- violating asymmetries in B'^ -^ ai(1260)^7r^ 
decays and determined the angle aeS [5]. In the absence of penguin contributions in these decay 
modes, this angle would coincide with the angle a of the unitary triangle of the Cabibbo-Kobayashi- 
Maskawa quark- mixing matrix [6]. Theoretical bounds on Aa = a — aeflf in these decay modes 
based on SU(3) flavor-symmetry have been derived in [7]. The rates of -B — > i^i(1270)7r and 
B -^ /Ci(1400)7r decays are experimental inputs to the calculation of these bounds. For the 
i^i(1400)'''7r~ decay mode |j there exists a published experimental upper limit at 90% confidence 
level (CL) of 1.1 x 10~^ [8]. Preliminary results for the branching fractions of the i^i(1270)"^7r~ and 
i^i(1400)^7r~ decay modes were obtained by the BABAR Collaboration on a sample of 384 million 
BB pairs [9]. In the following, we use Ki to indicate both i^i(1270) and i^i(1400) mesons. 

2 THE BABAR DETECTOR AND DATASET 

The results presented here are based on a sample of A'^^^ = 454.3 it 5.0 million BB pairs collected 
with the BABAR detector [lOj at the PEP-II e~^e~ asymmetric-energy storage rings. The e'^e~ 
center-of-mass energy y^ is equal to 10.58 GeV, corresponding to the T{4S) resonance. 

Momenta of charged particles are measured in a tracking system consisting of a silicon vertex 
tracker with five double-sided layers and a 40-layer drift chamber, both within the 1.5 T magnetic 
field of a solenoid. Identification of charged hadrons is provided by measurements of the energy 
loss in the tracking devices and by a ring-imaging Cherenkov detector. For lepton identification, we 
use the energy deposit in a CsI(Tl) electromagnetic calorimeter and the pattern of hits in resistive 
plate chambers (partially upgraded to limited streamer tubes for a subset of the data used in this 
analysis) intervalled with the passive material comprising the solenoid magnetic fiux return. 

3 ANALYSIS METHOD 

The B^ -^ K^TT^ candidates are identified from the K^ -^ K^n^n^ final state, with reconstructed 
mass mK-K-K in the [1.1, 1.8] GeV/c^ range. They are kinematically characterized by tit-es = [(•5/2 + 
Pj'-Pb)^/-E|- — p^]^ and /S.E = (£'r-E_B — pr'Ps — •s/2)/\/s, where [Eb, Pb) is the four-momentum 
of the B candidate, and [Ey^'Pt) is the e^e~ initial state four-momentum, both in the laboratory 
frame. We require mss > 5.25 GeV/c^ and |A_E| < 0.15 GeV. 

To reject the dominant e^e~ -^ quark-antiquark background, we use the thrust angle ^t be- 
tween the -B-candidate thrust axis and that of the rest of the event, calculated in the center-of-mass 
(CM) frame, and a Fisher discriminant J- [llj. The discriminant combines the polar angles of the 
S-momentum vector and the i?-candidate thrust axis with respect to the beam axis, and the zeroth 



Except as noted explicitly, we use a particle name to denote either member of a charge conjugate pair. 



and second moments of the energy flow around the i?-candidate thrust axis, calculated in the CM 
frame pjj. 

The resonant K^it^it~ system can receive contributions from several strange resonances in the 
selected range for mK-Kn, besides Ki mesons, such as K*{1410)^ {J^ = 1~), -fC*(1680)^ (1~)) and 
i^2 (l^SO)"*" (2"'"). Decays containing any of these resonances are characterized by different angular 
distributions. We define TC as the cosine of the angle between the direction of the primary pion 
from B decay and the normal to the plane defined by Ki daughters in Ki rest frame. We require 
l^l < 0.95 to reduce background from B^ -^ V^ir^ decay modes, where V^ is a vector meson 
decaying to K~^Tr~^Tr~. 

Background from B decays to final states with charm is suppressed by rejecting a signal candi- 
date if it has at least one track in common with a background B candidate, reconstructed in any of 
the B^ -^ D^TT^, B^ -^ D*^iT^, and B^ -^ D^-k^ background decay channels, with D meson mass 
within 0.07 GeV/c^ of the nominal value (if more than one such background candidates are recon- 
structed per event per background channel, the one with the highest B vertex fit x^ probability is 
chosen). To suppress background from B decays to final states with charmonium we calculate the 
invariant mass of the neutral vrvr combination of the primary pion from B decay with the opposite 
charge pion from Ki decay, and require that it is not consistent with any of the cc mesons J/ip, 
'ip{2S), rjc, ric{2S), Xco(l-P)) and Xci(l-P)- We also make particle identification requirements to 
identify pions and kaons, and veto muons, electrons and protons. 

The average number of candidates found per selected event in the data sample is 1.20. In case of 
events with multiple candidates, we select the candidate with the highest B vertex fit x^ probability. 
We classify the events according to the invariant masses of the K~^ir~ and TT~^7r~ systems in the 
K^ -^ K^TT^TT^ final state: events which satisfy the requirement 0.846 < mK-K < 0.946 GeV/c^ 
belong to class 1 {" K* region"); events not included in class 1 for which 0.500 < m^yr < 0.800 
GeV/c^ belong to class 2 ("p region"); all other events are rejected. 

A two-resonances, six channels i^-matrix model [U] is used to describe the resonant Ktttt 
system for the signal [13j. The notation is consistent with that used in [13]. The labels a and b in 
the following paragraphs refer to /fi(1400) and i^i(1270), respectively. The production amplitude 
for channel i = {{K* TT)s-wave, {K* TT)D-wave, P K,Kq vr, /o(1370) K,uj K} is given by 

F, = e'^- 5^(1 - iKp)r;p,, (1) 

j 

where 

T^ Jaijaj Jbijbj /n\ 

'^ ~ Ma-M + M,-M' ^ ' 

Si are offset phases ((5(_ft-*7r)g = 0), and P is the production vector 

p Jpalai . Jpbjbi /„\ 

' ~ Ma-M Mb-M' ^ ' 

The decay constants fai , fu and the X-matrix poles Ma and M\, are real. The elements of 
the diagonal phase space matrix p for the process Ki ^3 + 4, 3^5 + 6, where 4, 5 and 6 are 
long-lived pseudoscalar particles and 3 is a resonance, have been approximated with the form 



MM) = ^ 



m 1714 . ^ ^ . ^ , 

^(M -m* -mi + lA) 

m* + r?T,4 



1/2 

(4) 



where M is the mass of i^i, m,4 is the mass of 4, m* is the mean mass of 3 and A is the half width 
of 3 [13] • The parameters of K and the offset phases 5i are obtained from a fit to the intensity and 
the relative phases of the Kttt: channels, which were extracted by the ACCMOR Collaboration in 
a partial wave analysis of the data on the reaction K~p -^ K~tt~^tt~p accumulated by the WA3 
experiment [13]. For the fit to WA3 data we add a background term to the production vector |15j . 
The decay constants for the to K channel are fixed according to the quark model [13] . 

We express the complex production constants fpa and fph in terms of the production parameters 
C, = {0,4>)'- fpa = cos9, fph = sinOe^^, where 6 G [0,7r/2], and (j) G [0,27r]. In this parameterization, 
tan^ represents the magnitude of the production constant for the -ffi(1270) meson relative to that 
for the i^i(1400) meson, while (j) is the relative phase. 

Signal MC samples are generated by weighting the (/C+7r+7r^)7r^ population according to the 
amplitude X]i3^ajx(-^^'^^^ N)^*' where the term {K'^'K^'K^\i) consists of a factor describing the 
angular distribution of the (i^+vr^vr") system resulting from Ki decay, an amplitude for the reso- 
nant TT^vr" and K~^'k~ systems, and isospin factors, and is calculated using the formalism described 
in [16] . The branching fraction for Ki — > ojK is accounted for as a correction to the total selection 
efficiency. 

We use an unbinned, extended maximum-likelihood (ML) fit to extract the event yields Ug^r 
and the parameters of the probability density function (PDF) Vs,r- The subscript r = {1,2} 
corresponds to one of the event classes defined above. The index s represents six event categories 
used in our data model: 

• the signal S° -^ K^tt- (s = 1), 

• possible backgrounds from B^ -^ ai(1260) vr^ — > (vr 7r^7r^)7r^ (s = 2), 

• B° ^ L>-7r+ -^ {K+Tr-7r-)TT+ (s = 3), 

• B° ^ K*(1410)+7r- {s = 4), 

• B° ^ K*°7r+7r- + p^K+ir- {s = 5), 

• combinatorial background (.5 = 6). 

We perform a likelihood scan with respect to the parameters C, with 21 x 21 points. At each point, 
a simultaneous fit to the two event classes is performed. 

The signal and background PDFs are the products of the PDFs for independent variables. The 
signal PDFs for AE, ?7iESi and J^ are parameterized as the sum of Gaussian functions for the 
core of the distributions plus empirical functions accounting for the tails. The dependence on (^ of 
the selection efficiencies and the signal PDF for rnKmr are parameterized by means of templates 
modeled upon signal MC samples. Resonance production occurs in the non-signal B background 
and is taken into account in the PDFs. For the combinatorial background, we use polynomials, 
except for tties and J- distributions which are parameterized by an empirical phase-space function 
[TBJ and by Gaussian functions, respectively. The combinatorial background PDF is found to 
describe well both the dominant quark-antiquark background and the background from random 
combinations of B tracks. For all components, PDFs for TC are parameterized with polynomials. 

The likelihood £e for each candidate e belonging to class r is defined as Ce = ^g Ug^r 'Ps,ri^e', C; C); 
where the PDFs are formed using the set of observables Xg = {AE, niEs, ^, 'mKmr, T~(-} and the 
dependence on production parameters C is relevant only for the signal PDF. ^ represents all other 
PDF parameters. In the definition of Ce the yields of the signal category for the two classes 
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are expressed as a function of the signal branching fraction B as ni^i = B x N^^ x ei(C) and 
ni^2 = B X N^^ X e2(C)! where the total selection efficiency, eriC)^ includes the daughter branching 
fractions and the reconstruction efficiency obtained from MC simulation. 

The signal branching fraction is a free parameter in the fit. The yields for event categories s = 2 
and 3 are fixed to the values estimated from MC. The yields for the other background components 
are determined from the fit. The PDF parameters for combinatorial background are left free to 
vary in the fit while those for the other event categories are fixed to the values extracted from 
Monte Carlo (MC) simulation [17] and calibration B^ — > D~tt'^ decays. 

4 SYSTEMATIC STUDIES 

The main sources of systematic uncertainties are summarized in Table [TJ We repeat the fit by 
varying all the parameters in ^ which were not left fioating in the fit within their uncertainties, 
and obtain the associated systematic uncertainties. The signal PDF model excludes the fake com- 
binations originating from mis-reconstructed signal events. The biases due to the presence of fake 
combinations, or other imperfections in the signal PDF model are estimated with MC simulation. 
The finite resolution of the likelihood scan is also a source of bias. A systematic error is evalu- 
ated by varying the i^i(1270) and -ffi(1400) mass poles in the signal model, the parameterization 
of the intermediate resonances in Ki decay, and the offset phases 5i. Additional systematic un- 
certainty originates from potential peaking BB background, including B^ -^ ii'|(1430)^7r^ and 
B^ — > K^ {1680)'^ IT ~ , and is evaluated by introducing the corresponding components in the defi- 
nition of the likelihood and repeating the fit with their yields fixed to values estimated from the 
available experimental information [20]. We assign a systematic uncertainty due to yield variation 
in the B^ -^ ai(1260)^7r^ and B^ — > D~^ _ _7r~^ event categories. The above systematic uncer- 
tainties do not scale with event yield and are included in the calculation of the significance of the 
result. 

We estimate the systematic uncertainty due to the interference between the B^ — > K^tt~ and 
the B^ -^ K*^7r~^7:~ + p^K^iT~ decays using simulated samples in which the decay amplitudes are 
generated according to the results of this measurement. The overall phases and relative contribution 
for K*^-K^iT~ and p^K~^7r~ interfering states are assumed to be constant across the phase space 
and varied between zero and a maximum value using uniform prior distributions. We calculate 
the systematic uncertainty from the RMS variation of the average signal branching fraction and 
parameters. In the calculation of significance, this effect is assumed to scale with the square root of 
the signal branching fraction. The systematic uncertainties in efficiencies are dominated by those in 
track finding and particle identification. Other systematic effects arise from event-selection criteria, 
such as track multiplicity and thrust angle, and the number of B mesons. 

5 RESULTS 

Figure [1] shows the likelihood scan and the values of Bgg which minimize — In £ as a function of 

6 and (p. The absolute minimum occurs at 9 = 0.785 and (p = 0.942, and the signal branching 
fraction corresponding to that point of the scan is B{B^ -^ K^ir^) = (31.0 it 2.7) x 10~^. By 
interpolation between neighbouring points of the likelihood scan we extract 6 = 0.81 it 0.06 and 
(j) = 1.11 lb 0.28. The quoted errors on the branching fraction and production parameters C are 
only statistical and correspond to a 0.5 increase in — ln£. A second, local minimum is located at 
9 = 0.785 and (p = 3.454, and is associated to a 1.0 increase in — ln£. 
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Table 1: Estimates of systematic errors. For the branching fraction, some of these errors are additive 



(A) and given in units of 10 
combined in quadrature. 



others are multiphcative (M) and given in % . Contributions are 



Quantity 


B 


e 


^ 


PDF parameters (A) 


1.0 


0.01 


0.04 


MC/data correction (A) 


1.2 


0.05 


0.27 


ML Fit bias (A) 


0.6 


0.03 


0.02 


Scan (A) 


1.3 


0.04 


0.16 


Ki mass poles (A) 


2.2 


0.01 


0.36 


Ki offset phases (A) 


0.2 


0.01 


0.02 


Ki intermediate resonances (A) 


0.5 


0.00 


0.06 


Peaking BB bkg (A) 


0.8 


0.02 


0.27 


Fixed background yields (A) 


0.8 


0.04 


0.08 


Interference (A) 


5.9 


0.25 


0.52 


MC statistics (M) 


1.0 


— 


— 


Particle ID (M) 


2.9 


— 


— 


Track finding (M) 


1.0 


— 


— 


cos6It(M) 


1.0 


- 


- 


Track multip. (M) 


1.0 


- 


- 


Number BB (M) 


1.1 


- 


- 


Total (A) 


6.9 


0.26 


0.76 





Figure 1: Left: contour plot of — ln£ (no systematic effects included) in the 9,(f) plane; each line 
corresponds to a n^/2 increase in — ln£, with n = {1, 2, ...}, with respect to the minimum (indicated 
by a cross). Right: fitted value of Bgg, in units of 10~^ as a function of 6 and (j). 



A conservative estimate of significance is calculated from a likelihood ratio test A(— 21n£), 
assuming a x^ distribution with A^ = 3 degrees of freedom and minimizing the significance with 
respect to the production parameters {9, (p). Here A(— 21nvC) is the difference between the value of 
— 21n£ for zero signal and the value at its minimum for given values of Q (C represents the convo- 
lution of the likelihood with a Gaussian function representing additive systematic uncertainties on 
the branching fraction). We observe a non zero B^ -^ K^tt" branching fraction with significance 
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Figure 2: sPlot projections onto a) m-ES (class 1), b) rriES (class 2), c) AE (class 1), d) AE (class 
2) in the Kitt decay. Points represent on-resonance data, solid line is the signal fit function. 
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Figure 3: sPlot projection onto rnKmr for class 1 (left) and class 2 events (right). Points represent on- 
resonance data, solid line is the sum of the fit functions of the decay modes ii'i(1270)7r + iiri(1400)7r 
(dashed), i^*(1410)7r (dash-dotted), and i^*(892)7r-7r (dotted). Here the points are obtained without 
using any information about resonances in the fit, i.e. we use only rriES, AE, and J- variables, 
while for the normalization of the curves we use the signal yields obtained from the nominal fit. 



greater than 5.1 o". 

Figure [2] shows the distributions of AE' and m,ES for the signal events, obtained by the event- 
weighting technique (sPlot) described in [19]. For each event, a weight to be signal or background 
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Figure 4: 68 % (dark shaded zone) and 95 % (light shaded zone) probabihty regions for 9 and 
(top), 9 (bottom-left) and (f) (bottom-right). 



is derived according to the results of the fit to all variables and the probability distributions in the 
restricted set of variables, in which the projection variable is omitted. Using these weights, the 
data is then plotted in the projection variable. We show in Figure [3] the projection onto mKmr- 

The experimental two-dimensional likelihood C for 9 and (j) is convoluted with a two-dimensional 
Gaussian that accounts for the systematic uncertainties. In Figure [4] we show the distributions we 
obtain for 9, cj) and 9 vs. (j) (the 68% and 95% probability regions are shown in dark and light shading 
respectively, and are defined as the regions which satisfy C(r) > Cmin and (r/^\~^r C{r)dr = 
68% (95%), where r is the projected set of variables). The condition £(r) > Cjj, where £// is the 
value of the likelihood evaluated at the position of the second, local minimum in Figure [H defines 
a 48% probability region, with systematic uncertainties included, on the 9 vs. (p plane. 



6 CONCLUSIONS 

We measure the branching fraction 



BiB° 



K+TT- 



) = (31.0 ±2.7 ±6.9) X 10 



-6 



with significance greater than 5.1 a. The first error quoted is statistical and the second systematic. 
The value of the branching fraction measured in this analysis is consistent with preliminary results 
obtained by BABAR Collaboration ^, and is to be compared with the naive factorization [21 [3] and 
QCD factorization [4j estimates, of order 10^^. 
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For the production parameters we obtain 

0.25 <e < 1.32 
-0.51 <(j)< 4.51 

at 95% probability. This analysis represents the first attempt to measure the relative phase between 
the production amphtudes of i^i(1270) and i^i(1400) mesons in B decays. 
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